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ABSTRACT 



Propagation of extremely short pulses of electromagnetic field (electro- 
magnetic spikes) is considered in the framework of the total Maxwell-Duffing 
model where anharmonic oscillators with cubic nonlinearities (Duffing model) 
represent the material medium and wave propagation is governed by the 1-d 
bidirectional Maxwell equations. This system of equations has a one pa- 
rameter family of exact analytical solutions representing an electromagnetic 

spike propagating on a zero or a nonzero background. We find that the total 
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Maxwell-Duffing equations can be written as a system in bilinear form and 
that the one-soliton solution of this system coincides with the steady state 
solution obtained previously. 
PACS: 42.65. Tg 
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1 Introduction 

Ultra short nonlinear pulses of an electromagnetic field, which contain as few 
as one half optical cycle, have recently attracted a great deal of attention 
PP-[H]. The description of the evolution of the electromagnetic field was 
based on Maxwell equations or the subsequent wave equation. To describe 
a medium where electromagnetic waves propagate, one frequently uses an 
ensemble of oscillators or resonant atoms. If the oscillators are linear we 
obtain the important Lorentz model, which has been very useful to describe 
the propagation of an electromagnetic wave in a linear medium. The simplest 
generalization of the Lorentz model is obtained by adding an anharmonic 
term to the equation of the oscillator. This leads to the Duffing model in the 
case of a cubic nonlinearity [HI El IE] ■ In the framework of the Duffing model 
it was shown that solitary pulses of a unidirectional electromagnetic field 
("electromagnetic bubbles") having an amplitude comparable to the atomic 
field and a duration down to ~ 10~^^ s may be expected. These bubbles 
propagate without dispersion as stable solitary waves. Very short pulses of 
this kind will be referred to as extremely short pulses (ESP). 
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In all these works, the oscillator represents the response of the high- 
frequency electronic degree of freedom to the electromagnetic field. The 
steady state ESP, which have been found heretofore, can be represented by a 
moving wave packet with a zero dominant frequency. Some authors call such 
pulses video pulses El to emphasize the difference from the second kind 
of ESP with a non-zero dominant frequency. These (optical) ESP are con- 
sidered in Ref. j^. Hereafter we will consider ESP's of the first kind. Since 
the ESP spectrum is concentrated at low frequencies, ion oscillations may 
also give a considerable contribution to the full polarization of the medium. 
Nevertheless, we will neglect the ion component in the material response (the 
propagation of femtosecond pulses in a medium with a nonlinearity deter- 
mined by both electronic and ionic (Raman-scattering) degrees of freedom 
was considered in Ref. pHl). 

As the Duffing model is the simplest and most fruitful generalization of 
the Lorentz model, its investigation is very attractive. An objective of the 
present work is to study the propagation of linearly polarized ESPs in a non- 
linear dispersive medium modeled by an anharmonic oscillator characterized 
by the cubic nonlinearities. The evolution of the electric field of ESP will be 
considered on the base of the Maxwell equations without any unidirectional 
reduction. 

The paper is structured as follows. The model is derived in section 2 
together with dynamical invariants or integrals of motion. Two families of 
moving ESP analytical solutions are given in section 3. In section 4 the 
bilinear form of these equations will be presented and we conclude in section 
5. 
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2 Constitutive model 



To consider the propagation of an extremely short electromagnetic pulse 
propagation in a nonlinear dispersive medium we need to use the total 
Maxwell wave equation and some model for the medium. Let us consider 
the standard Lagrangian density for an electromagnetic field taking into ac- 
count the interaction with an ensemble of nonlinear oscillators: 

^ = ^^^1 - + A'nY. {^^^a, - \niXl - \K,^Xt - ^AX,,} (1) 

Here we consider a plane electromagnetic wave, propagating along the 
2-axis and represented by the vector potential A. We use an anharmonic 
oscillators model (Duffing model) to reproduce the electronic response of an 
atom located at the spatial point indicated by the symbol a. The electrons are 
considered as particles in a potential well characterized by the displacements 
from their equilibrium position Xa- They oscillate with their own frequencies 
uja and are influenced by the electromagnetic field. In expression ((H) is 
the electric charge of the electron and n^a are anharmonicity coefficients. 
Hereafter, we will use m as a symbol for this effective mass, which accounts 
for the local Lorentz field effect. The partial derivatives we will denote as 
df/dt = f^t, df/dz = and so on. 

The application of the variational procedure to the action related with 
the Lagrangian density ((H) yields equations 

A^zz - c'^A^tt = {4:7ie^/co)'^Xa,t, 

a 

mXa,tt + mulXa + Ksa^a = (e*/co)A,i 
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If one introduces the strength of the electric field E = Cq ^A^t, then the 
constituent equations of the model under consideration can be written as 

- c-^E^u = {^T^/cl)P,tt (2) 
Xa,tt + ^IXa + {>^za/m)Xl = {eJm)E (3) 

The polarization P of the nonlinear medium is P = X^a e*Xa . 

Let us consider the case of a homogeneous broadening medium where all 
atoms have the same parameters, in particular cj^ = cuq- Then we can write 
the polarization as P = n^e^X, where ha is the density of the oscillators 
(atoms), and the index of the atom can be omitted. 

We rescale the variables and fields as C = zujq/cq, t = uj^t and e = E/Eq, 
q = X/Xq, where 

Eo = mujlXo/e, = mc^||e;^(2/i/|/^3|)'/^ Xq = {2fimujl/\K3\f\ 

and we will also use the following parameters 7 = oUp/uo, 2fi = K^XQ/mul, 
where Up = {ATmAel/mY^'^ is the plasma frequency. In terms of the rescaled 
variables equations Q and Q take the form 

e,cc - (^,TT = iq,TT , q,TT + q + '^m^ = le. (4) 

These equations will be named the total Maxwell- Duffing equations (or TMD- 
equations). To this system should be added the initial conditions: 

e(C = 0,r) =eo(r), e,.(C = 0, r) = e[,(r), 

and the boundary conditions 

qiC,T) = qAC^t) = 0, at r ^ ±00. 
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related to the evolution of an initial electromagnetic solitary wave in a non- 
linear dispersive medium. 

Note that the system of equations Q can be rewritten in alternative form 
by introducing the new auxiliary field variable b by the relation de/d( = 
db/dr. In this case the TMD-equations take the following form: 

e,c = &,r , fe,c = e,r + 1P , q,T=P , P,T + q + = 'je. (5) 

It should be noted that the system of equations (jlj) can be derived as the 
Euler-Lagrange equations from the action functional 



S = J ^[q^ a]dTdC, 



where now the new Lagrangian density is 

1 (daV 1 (daV 1 ( dqV 1 ^ u ^ da 

Applying the variational procedure to the action 5* yields equations 
d'^a d'^a dq d'^q o da 

Identifying potential for the field e, so that e = a^r, makes these 

equations identical to Eqs. (jlj). 

From the Lagrangian density (jHj) we can get the density of moments of 
the fields a and q: 

7r.(C,r) = |^ = ac(C,r) = 6(C,r) , vr,(C, r) = |^ = 0. (8) 

The second expression in (jHj) indicates that we have a degenerate Lagrangian, 
which leads to a constrained Hamiltonian system {TTqlC, r) = is a primary 
constraint) [T2JIISj- 
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Now one can get a canonical Hamiltonian density using the Legendre 
transform 

dC 1 

'H = ^ = 2 ('^^'^ ^'^^ ~^ ~^ ~ ^^"^ ^ 27ga,^) . (9) 

The integration of this expression with respect to r leads to the canonical 
Hamiltonian which is an integral of motion 

1 °° 

H = ^ J (e^ + h'' + ql-q^ - fiq"^ + 2-fqe) dr. (10) 

— oo 

It is worth noting that there are two additional integrals of motion, which 
follow from the TMD-equations (0): 

oo oo 

h= J e{C,T)dT= J a^(C,r)dr = a(C,r = oo) - a(C,r = -oo), (11) 

— oo — oo 

and 

oo oo 

l2= J b{C,r)dr= J n,{C,r)dT (12) 

— oo — oo 

The magnitude of the first integral is defined by the boundary conditions only 
so that it can be interpreted as a topological charge in the Maxwell-Duffing 
model. The second integral is the canonical moment in this model. 



3 Steady state solutions 

Let us look for solutions of the TMD-equations as traveling waves with a 
non-varying profile. First consider solitary waves on a zero background, i.e. 
assume the following boundary conditions at r — > ±oo: 

e(C,r) =e,.(C,r) = , g(C, r) = g,.(C, r) = 0, 
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Substituting e(C, t) = e{T — (/v), q{Cj'^) = into equation (jD) and 

taking into account the boundary conditions one finds that the first equation 
of (jH) results in e = a;g/7, where a = 7^f^/(l — v^). The second equation 
of (jU can be transformed into the following ordinary differential equation in 
the variable T = r — C,/v 

^-(a-l)g + 2/xg3 = o. 

A non-singular solution of this equation exists only if a > 1 and > 0. In 
this case integrating results in the following expression |2l El El 

<lst{C,,T) = ±^{a- l)//isech \/(« - 1) (r - C/v - Tq) . (13) 

Here the integration constants tq and a are the parameters of this steady 
state solitary wave. The strength of the electric field of the ESP is given by 
the following formula 

estiC, r) = ±a7~V("-l)/'"sech y {a - 1) (r - (/v - To)] (14) 

This solution describes the propagation of an electromagnetic spike with 
positive (+ sign ) or negative (— sign ) polarity. The condition for existence 
of this solution leads to the limitation of the velocity: 1 < i; < (1 + 7^)^^/^. 

Now let us consider solitary waves on a non-zero background, i.e. we 
assume the following boundary conditions at r — > ±00:: 

e(C,^)=eo, qiC,T) = qo , e,^(C,r) = , g,^(C,r) = 0, (15) 

7eo = qo + (16) 
Introduce the following fields u{(, r) = e{(, r) — cq and f{(, r) = q{(, r) — go- 
Instead of equations (jl} we obtain the new system of equations 

^-^ = 7^ , ^ + (l + 6/igo)/ + 6/^W +2/./ =-i—. (17) 
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Under the assumption that u{(,t) = u{t — fiC^'^) = fij ~ Ch) the 
first equation leads to the relation e = aqj^^ and the second equation takes 
the form 

-± = {a-l- Q^iqDf - Gfigof - 2fif. (18) 
From (fTH|) and the boundary conditions (fT^ it follows that 

= (« - 1 - 6/xgo')f - ^mf - l^f- (19) 

A non-singular solution of this equation exists when ai = a — Gfiq^ > 1 and 
/X > 0. This condition leads to the limitation of the velocity of the ESP 

2 2 1 + 

1 <v'' <vi - ^ 



1 + 7^ + 6/igo ' 
The solution of equation can be written as 

±(ai - 1) 



^(2/igo)2 + /i(«i-l)2cosh[(«i - l)i/2(r - (/t; - tq)] ± 2/igo 

(20) 



The steady state pulse of electromagnetic wave propagating on a constant 
electric background is then given by 

y(2/xgo)' + /i(ai - 1)' cosh[(ai - 1) V2(^ _ (^/y _ r^)] ± 2/igo 

(21) 

Unlike the case of ESP propagation without background, here ESPs of dif- 
ferent polarities have different amplitudes. 



4 Bilinear form of the total Maxwell- 
Duffing equations 

If the following substitutions 



(22) 
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are used, then equations (0) can be rewritten as 

^^D^{g-h)~-^^DAc-h)=0, 

^D^{c ■ h) - ^D^g ■ h) = l^Drif ■ h), (23) 

where the Hirota /^-operators i5^(a ■ b) = a^cp — ah^c_, Dr{a ■ b) = a^rb — ab^r 
have been introduced [H]. To derive the last equation we follow the rule 



-Diif.h)-{-D^Ah-h). 



Multiplying the first equation by /i^, the second equation one by h^, the third 
one by h"^ and collecting the terms of order h"^, we can write the resulting 
equations as a system of bilinear ones 

D^{g-h) = D^c-h), 

D^{c-h) = D,{g-h)+jD,{f-h), (24) 

DUf-h) = {^g-f)h, 

Dl{h-h) = 2fif 
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We will use the usual method [HI HH] to solve equations (^1]) by writing 

g = gie\ c = cie\ f = he\ h = I + h^e' + h^e^', (25) 

where 6 = k( — Qt. Substituting this into results in the equations: 

kgie^ - kh2gie^^ + flcie^ - Vth2Cie^^ = 0, 
kcie^ - khiCie^'^ + O^ic^ - nh2gie^^ + ^7/16^ - n^h^he^^ = 0, 
n^f,e' + n'hf^e''-{Jg^-f,){e' + he'' + h2e'') = 0, 

Equating the coefficients of the different powers of to zero, one obtains the 
system of equations that define Ci, fi,gi, hi h2 : 

kgi + Vtci = , kci + n{gi+'jfi) = 0, 
n^fi-hgi-fi) = , 4fi2/,2=^/2 , h, = 0. 

From these relations one can get 

h = , /i2 = /i/i'(2fi)-2, 

ci = -r\k/n){i + n^)h , g^ = r\^ + ^')fl, 

and the "dispersion relation" 

It should be pointed out that in the low frequency limit this formula yields 
/c^ = + 7^), while in the high-frequency limit it yields k'^ = Q"^. 
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Thus, we found a one-soliton solution of the bihnear equations (j24p which 
can be written as 

g = 7-^(1 + ^^')/le^ 

c = -^-\k/n){i + n'')he\ 

f = he' , h = l + fif,i2n)-'e''. 

These relations yield a solution of the TMD-equations (jH)) which is consistent 
with the steady state solution obtained earlier: 

^ " l + /i(//2^])v«' ^ ' 



he' 



l + ^(//2^])2e2«■ 

If we introduce a new constant of integration exp6'o = ±yU^/^/i/2n, then the 
one-soliton expression for the normalized electric field of ESP can be written 

as 

esoiiC, r) = ±i 29 
\ fi^/'^'y J cosh(t^ + t^o) 

The velocity of this soliton is defined as v = Q/k. Taking into account the 
dispersion relation (|26|l one can obtain 



1 + 



The magnitude of this velocity lies in the interval (1, (l+7^)~^/^). This is the 
same result as for a steady state pulse obtained above. Using the expression 
for velocity we can see that a = 1 + Q"^, and the expression for the steady 
state pulse (fT^ coincides with expression ((221) • 
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5 Conclusion 



We have introduced and analyzed a model for the propagation of extremely 
short unipolar pulses of an electromagnetic field in a medium represented by 
anharmonic oscillators with a cubic nonlinearity. The model under considera- 
tion takes into account the dispersion properties of both linear and nonlinear 
responses of the medium. This model is the simplest generalization of the well 
known Lorentz model used to describe linear optical properties in condensed 
matter. The cubic nonlinearity is the first type of anharmonic correction to 
the Lorentz model and it results in the Duffing oscillator. Here we consider 
the total Maxwell-Duffing model in detail. The Lagrangian picture of the 
TMD model was considered and three integrals of motion were found. Two 
families of exact analytical solutions, with positive and negative polarities, 
have been found as moving solitary pulses. The first kind of steady state 
ESP is an electromagnetic spike propagating in a nonlinear medium. It was 
discussed early in 13 El CHI- A new kind of steady state ESP is an elec- 
tromagnetic spike propagating on a non-zero electric background. There are 
both bright and dark ESP. Unlike the ESP on a zero background, here pulses 
of different polarities have different amplitudes. 

We found that the TMD equations can be represented in bilinear Hirota's 
form. In the case of a zero background the one-soliton solution of the bilinear 
equations was obtained. It coincides with the expression of a steady state 
ESP. There are many examples when the bilinear form of the nonlinear evolu- 
tion equations leads to the existence of two-soliton solutions without ensuring 
complete integrability or/and the existence of A^-soliton solutions. However 
in this particular case we were not able to obtain a two-soliton solution of 
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the TMD equations. 
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